Theory of a double-dot charge detector 



00 

o 
o 

(N 
X) 

(N 



Tamas Geszti and Jozsef Zsolt Bernad 
Department of the Physics of Complex Systems, 
Eotvos University; H-1117 Budapest, Hungary 
e-mail: geszti@complex.elte.hu 

A double quantum dot charge detector, with one dot Coulomb coupled to the electron to be 
detected and the other modulated by a time-dependent plunger voltage, is analyzed in a minimal 
model. The signal and noise of the detector are calculated by a standard master equation and 
MacDonald formula technique. We find a dip in the noise spectrum at the double Rabi frequency 
of the double dot, defining the bandwidth available for detecting charges in motion. 
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Mesoscopically narrow passages of electrons, control- 
ling the flow of electric current by potential barriers just 
below or above the Fermi level, are obviously very sen- 
sitive to external electric fields, which suggests their use 
as charge detectors. Most often, a mesoscopic trap or 
"island" carrying the charge is Coulomb coupled to a 
quantum point contacti or a quantum dot - in that ar- 
rangement called a single-electron transistor (SET),^, to 
offer a charge sensitivity sufficient to detect individual 
electrons. Similar tools are being explored as readout 
devices for charge qubits^ with some results emerging 
about the quantification of the amount of information 
accessible in the readout process 4 

Most recently, double quantum dots (DQDs))^ here 
also called double island single-electron transistors, have 
been used for charge detection,^ the main advantage be- 
ing their relative immunity against noise of various ori- 
gins. The operation is based on the highly coherent inter- 
dot tunneling, manifest through Rabi oscillations,^ pro- 
viding sharp boundaries on the current-voltage charac- 
teristic of DQDs.^. This new kind of device has been 
theoretically analyzed by Tanamoto and Huj^ through a 
model mainly focusing on the complications from weak 
Coulomb blocade, and specific aspects related to quan- 
tum information, also discussed in Ref. ITol . 

In most of the experimental work cited above, the de- 
vice is operated at radio frequency, drawing on analogous 
work on radio-frequency SETs,^^ efficiently avoiding the 
band of low frequencies, strongly contaminated by 1// 
noise, which is ubiquitous, and hard to analyze theoret- 
ically. On the contrary, shot noise, extending to much 
higher frequencies, is open to analysis by standard tools, 
which is also our main concern here. 

In the present paper we discuss charge detection in 
DQDs, supplemented by the possibility of fast time con- 
trol, as sketched in Fig. [1] The island trapping the elec- 
tron to be detected would be placed next to the first 
compartment of the DQD, so the electron would detune 
the tunneling resonance set by the two bias plungers:^^ 
Time control, as inspired by the experiment of Nakamura 
et alr^ could be achieved by a synchronized modulation 
of the pump voltage sending the electron onto the is- 
land, and the second plunger of the DQD, driving the 
device through a work point. A variety of delay times 
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FIG. 1; Double quantum dot (DQD) electron detector. The 
electron to be detected is driven by pump voltage Vpump onto 
the island (trap) D* , then back; its presence or absence is 
reflected in the current through the DQD. The island is 
Coulomb coupled to compartment Di of the DQD, set to a 
fixed plunger voltage Vp^ ; plunger voltage of compartment 
D2 is used to time-gate the detector. 



and plunger pulse shapes would allow a detailed analysis 
of the dynamics of one-electron detection. 

In calculating signal and noise of the DQD detector, 
we evaluate the reduced density matrix of the coupled 
DQD-|-trap system^^i by means of a standard Markovian 
master equation approach (see, e.g., Refs ITsllTGllTTh . sup- 
plemented by the possibility of resolving the dynamics 
according to the number of achieved tunneling events;— 
The latter being slow, their fast screening in the exter- 
nal metallic circuits can be approximately treated by 
the Ramo-Shockley theorem in its symmetric formJ^ 
Frequency-dependent stationary noise is evaluated by 
means of MacDonald's formulajAS which is a shortcut 
avoiding the more familiar Wiener-Khinchine analysis of 
time correlation functions. 

The physical time and energy scale of the processes en- 
visaged is set by the Rabi frequency. For a typical DQD 
formed by metal contacts on top of a GaAs/AlGaAs 
heterostructure^ it is of the order « 10^° Hz. To 
preserve coherence, tunneling rates should be set by gate 
voltages to the same order of magnitude, thereby assuring 
time control on the ns scale, easy to follow by external 
electronics. 

Coulomb blockade allows one to consider the low- 
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temperature dynamics of the DQD + trap system on 
a truncated basis, stretched by states of or 1 extra 
electron above a neutral background on each dot and 
the trap. A further step of truncation is brought about 
by the interdot Coulomb repulsion which for the geom- 
etry of the experiment of Hayashi et al— amounts to 
« lO'^ h ft, effectively excluding states for which both 
dots are occupiedJiii? 

The above considerations justify the use of a minimal 
model in which the retained basis vectors are denoted as 
\a) = 1000) , \b) = |1GG) , |c) = 1010) , \d) = |0G1) , |e) = 
101), I/) = |011); the three occupation numbers be- 
longing to left dot, right dot, and trap respectively. In 
addition, each of those three parts is coupled to a sepa- 
rate metallic contact, acting as a fermion reservoir. ^'^ The 
dynamics of the whole is generated by the Hamiltonian 
H = Hdqd + Ht + Hint, with 

HoQa/fi eia|ai + £202*^2 + ^{a\a2 + ajai) 
+ wibjbi + Wrblbr + Y K^lai + h.c. 

I r I 

+ YKbla2 + h.c, (la) 

r 

Hrp/h — ^ Wpdpdp + erc^c + ^ A*dpC + h.c, (lb) 

p p 
Hint/h = vaXaic^c (Ic) 

Reduced dynamics of the DQD + trap subsystem are ob- 
tained by carrying out thermal averaging over the states 
of the contacts in the initial state. In the above formu- 
las, hti {i = 1,2) is the energy of an electron occupying 
the single state allowed by Coulomb blockade in the ith 
quantum dot, hfl is the amplitude of tunneling between 
the two dots, hXi and hX,, are the respective tunneling 
amplitudes from left reservoir to left dot and from right 
dot to right reservoir; hv is the Coulomb interaction ma- 
trix element characterizing the coupling of the trap to the 
first dot, hXp is the tunneling amplitude between pump 
and trap, and her is the energy level of the trap; hwi and 
hwr are the one-electron energies in the left and right con- 
tacts, and hwp in the pump, is the one-dot electron 
annihilation operator, bi and br are those for the respec- 
tive contacts, c is that of the electron trap, and dp that 
in the pth one-electron state of the pump. 

The one-dot levels Ci can be modulated in time through 
plunger voltages Vp. to scan near-resonance conditions. 
Introducing the notation (5 ei — £2, 
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is the condition of resonance, assuring that if the trap 
is occupied by one electron, the DQD detector exhibits 
maximum dc transmission, i.e., maximum signal with re- 
spect to the empty trap case. 

As mentioned above, in conjunction with a Ramo- 
Shockley treatment of screening, it is convenient to de- 
compose the density matrix according to the total num- 



ber A'' of tunneling events occurring at both external con- 
tacts of the detector, which immediately furnishes the 
Ramo-Shockley screened current asi^ 
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where 



with the notations 

— \Paa J Pbb ' Pbc ' Pcb ^ Pec ) J 
En — ^Pdd ' Pee i Pef ' Pfe ' Pff ) ' 

and 

v=il, 1,0, 0, If. 
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The vectors x^^it) and Vj^it) are determined by the solu- 
tion of a system of master equations, with the appropri- 
ate initial conditions, which are by no means trivial.— 

We restrict ourselves to the zero-temperature case; 
then the Fermi levels of the contacts, which control the 
detector and the pumping of the trap, respectively, do 
not appear explicitly in the calculation, apart from distin- 
guishing pumping-in and pumping-out periods, according 
to whether the pump electrode Fermi level is above or 
below the trap level hex- Starting with Hamiltonian ([1]), 
and following any of the equivalent standa rd pro cedures 
used, e.g., in the papers listed under Refs. [1111113, one 
invariably arrives at a system of Markovian master equa- 
tions for the components of the density matrix, which 
can be written in the form 



= A: 
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for pumping-in, 
for pumping-out. 



(7a) 
(7b) 
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and we have introduced the matrices 
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Here and Tn are the familiar Fermi golden rule 
rates of dissipative tunneling transitions from left con- 
tact to left dot and from right dot to right contact, 
respectively]^ F^ is the rate of transitions from pump- 
ing contact to trap and vice versa. Damping and dcphas- 
ing of coherent processes in the DQD detector are con- 
trolled by these tunneling rates, microscopically rather 
ill defined because of their high sensitivity to defects and 
impurities, but - as already mentioned - easy to control 
experimentally by gate voltages. 

Trivially, the same system of master equations is 
obeyed by the vectors x, y composed of the unresolved 

density matrix elements pij = p\J^ . It is easy to find 
their steady-state solution pfj""^; then through Eq. ([3|) 
we get the steady-state current^^ 
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FIG. 2: Current signal of the detector as function of time, 
dimensionless scales being set by interdot tunneling ma- 
trix element Q and electron charge e. Parameters T/^l — 
2/3, Fr/n = 1 are tunable by plunger voltages. For coupling 
of relative strength v/Q = 10/3 [moderately strong coupling; 
see Eq. (|12|l ] which is fixed by construction, the DQD is bi- 
ased to S — —11/3 to achieve high charge sensitivity. Charge 
is pumped onto the island at i}t = 20 and pumped out at 
Qt = 32; the detector is shut down at fit = 50 by shifting 
bias S to -8. 



where 



6 for empty trap, 

S + ly for filled trap; 



(11) 



in addition, we have introduced the asymmetry param- 
eter a = (Fi — Tii)/{Tl + Ffl) and the mean tunneling 
rate F = (Fl -t- Ffl)/2. 

The zero-frequency signal of the detector is defined by 
the difference between the current with filled and empty 
trap, the detector being active all the time. Large asym- 
metries (i.e., a close to unity) cut down the current, 
which is the signal of our detector. Although all the 
subsequent analysis can be carried through for arbitrary 
values of a, the gain thereby is negligible, therefore in 
what follows, we restrict ourselves to the symmetric case 
a = 0. 

Noting that (5 is a "soft" experimental parameter, easy 
to adjust as desired, we observe that the sensitivity of the 
current to the presence or absence of the coupling term 
V essentially depends on whether that coupling is strong 
(j/ ^ A) or weak {ly <^ A) with respect to 



(12) 



The conclusion for the experimenter is this: to achieve 
maximum signal, strive at strong coupling in the above 
sense, and set the bias to S — --v. Taking again the 
Hayashi et al. experiment^ for reference, in all practical 
cases we are in the strong-coupling limit; we note that 
the actual Coulomb coupling strongly depends on the 
geometry and can be modulated by shifting the island. 

Detector performance depends on noise as well; in 
particular, on shot noise, caused by the discreteness of 
electron transport under the joint control of tunneling 
and Coulomb blockade.— Time-dependent noise analy- 
sis, with optimized filtering of short and noisy time- 
gated pulses responding to time-dependent charging of 



the trap, is a difficult issue of signal processing. For the 
present paper, we neglect all those complications implicit 
in Eqs. ([7]), and restrict ourselves to the case of constant 
charge on the trap. Then it is straightforward to calculate 
the frequency-dependent steady-state noise spectrum of 
the DQD detector, using the now standar d^^'^^ technique 
based on MacDonald's formula'^^'-P 



sin{LUT). (13) 



Here Q{t) is the total charge carried until time t, deter- 
mined through the time integral of Eq. Q , which gives 



dt 
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in which the long-time linear time dependence of the first 
term on the right-hand side is canceled by the last term. 
Then the desired noise spectrum is evaluated in the form 



S(lu) = Lu— \M(^iuj) - M(iuj)] , (15) 
4i 



M{z)= dt e-^^^N'^pNit), (16) 



nih. 



to which, as expected, the stationary-current term of 
Eq. gives no contributioni^ Accordingly, we pro- 

ceed through the Laplace transform solution of Eqs. ([7]), 
using the initial conditions^* 
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Having taken the Laplace transform, Eqs. ([7]) furnish 
an iteration in N from to cxd, presenting the Laplace 
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transform of Eq. p4|) as a matrix series, which can be 
summed up in a closed form to give 



M{z) 



K + K^ 



(zI-A)-^ xo(0) (18) 



with K — (zl — A) B, the parameter S in matrix A 
[see Eq. (j9])] being r eplaced by S [Eq. ^] to cover 
both cases of empty and filled traps. The matrix inverses 
can be evaluated analytically, to finally give the scaled 
formula 



^(wir,^,^) = e^n s 




(19) 



where s{x \y, z) — u(x, y, z)/v{x, y, z), with 

u{x,y,z) ^Ay (l6a;^ + 8x^ (7?/^ -4 (4 + 2^)) 

(^572/-* + 16 (4 + 2^)2 _ (46 + 1122^ 

+2?/4 (y^ + S?/^ {-l + z^) +16 {b + Qz^ + z^)) 
Wy^ (l9/-4t/2 (37+10 2^) 
+ 16 (44 + 23 2^ + 3 2''))) , 



(20) 



and 



v{x,y,z) = 16 2:^ + y'' ((y^+4 (3 + 2^))) +8a;^ {5y^ 



x'^ (33y* + 16(4 + 22) -sy- 



„8 , y4 //„,2 , Wq I .,2 
-4 (4 + 2^)) 

(32 + 7 22)) +2a;2?/2 (Sy-^ + y^ (20 - 8 2^) 
+ 16 (20 + 9 2^ + 2^)) . 

(21) 

The results^ are displayed in Fig. [31 The most con- 
spicuous feature is a dip appearing at frequencies uj 
slightly below 2fl, the frequency of undamped Rabi os- 
cillations. The dip is visible under near-resonance condi- 
tions (5 ~ 0, corresponding to peak signal. We think the 
qualitative reason is this. Shot noise is generated when 
an electron enters or leaves the DQD at an external con- 
tact, enforcing partition of the current. One oscillation 
period after entering from the left, the electron is just 
back to the left dot, so it cannot leave right. That kind 
of history is quantified by the solution of the system of 
master equations ([7]). 

In conclusion, our "minimal model" offers some in- 
sight into the potentialities of charge detection using a 
double quantum dot. That refers particularly to the 
noise characteristics: 1/f noise being excluded by radio- 
frequency operation, intensive shot noise is expected to 
appear above the Rabi frequency 257, and if not cut by 
the RC time in the external electronics, it can be effi- 
ciently reduced by low-pass filtering. That sets the in- 
trinsic speed of detection to the time scale f^"^, which in 



turn depends on geometry and gating voltages; according 
to Ref. 7, the DQD detector itself can be good down to 
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FIG. 3: Scaled noise spectrum of the double-dot charge de- 
tector (see text), as a function of relative frequency lu/Q and 
relative bias 5/Q, [see Eq. with relative tunneling rate 

fixed to r/n = 1.2 



the nanosecond scale. The practical limitation in reach- 
ing that performance is typically the external circuit. 

It is worth mentioning that such arrangements might 
serve as tools for studying fundamental issues re- 
lated to the dynamics of the quantum measurement 
processi^i^^i^ Indeed, measuring two-detector correla- 
tions between detectors differently gated in time may 
give a new chance to experimentally approach the "col- 
lapse of the wave vector," which has so far notoriously 
resisted being resolved as a physical process with its own 
dynamics. Although a disadvantage in the qubit read- 
out issue, in this context it may turn into an advan- 
tage that the operation of fast single-electron detectors 
is still slower than that of their single-photon counter- 
parts, therefore there is more real chance to achieve the 
necessary time resolution for the purpose. Deterministic 
single-electron sources, based on various pump or turn- 
stile constructions^^ have existed for some time, offering 
the possibility to time-lock detectors to them. What one 
may happen to see is nonstandard precursor fluctuations 
in the detector current, as the random choice inherent in 
the quantum measurement process develops in time. 
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